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STATISTICAL  ANALYSIS  OF  DETECTION  PERFORMANCE 
FOR  LARGE  DISTRIBUTED  SENSOR  SYSTEMS 


1.  INTRODUCTION 


This  report  develops  analytical  expressions  for  the  detection  performance  of  a  sensor  system 
composed  of  a  large  field  of  distributed,  short-range  detectors.  If  a  field  can  be  deployed  slowly 
and  sensors  anchored  to  a  specific  location,  the  performance  can  be  analyzed  using  traditional 
search  theory  for  a  large  number  of  sensors.  Furthermore,  if  the  number  of  sensors  is  large 
enough  so  that  the  field  coverage  is  very  dense,  visual  examination  of  the  resulting  distribution  of 
“contacts”  will  remove  false  alarms  and  detection  performance  of  the  field  can  be  analyzed  by 
examining  the  detection  performance  of  an  individual  sensor.  Neither  one  of  these  analysis 
techniques  is  appropriate  if  the  field  is  sparsely  populated  in  a  random  manner,  and  that  is  the  type 
of  sensor  field  analyzed  in  this  report. 

The  analysis  in  this  report  examines  the  issues  of  detection  performance  of  large  distributed 
systems  of  proximity  sensors.  The  effects  of  sensor  field  distribution  (both  uniform  and 
non-uniform),  sensor  performance  characteristics,  and  basic  tracking  strategy  on  the  performance 
of  the  overall  field  are  examined.  The  performance  is  analyzed  in  terms  of  overall  probability  of 
detection  and  probability  of  false  alarm  during  some  operation  time.  The  report  begins  with  a 
discussion  of  sensor  performance,  incorporating  range  dependency  in  ideal  and  non-ideal 
scenarios.  The  next  section  discusses  field  detection  performance  for  both  stationary  and  moving 
targets.  The  third  section  of  the  report  discusses  some  field  design  strategies  based  on  the 
analysis.  The  report  concludes  with  a  summary  of  the  primary  analysis  results. 


2.  DISTRIBUTED  AREA  COVERAGE  MODEL 


The  modeling  of  the  detection  performance  of  a  distributed  sensor  system  requires  three 
uncertain  quantities  that  are  modeled  with  random  variables:  individual  sensor  performance, 
sensor  distribution  in  the  field,  and  target  location/movement.  These  parameters  are  treated  as 
random  with  some  known  characteristics.  For  the  purposes  of  this  analysis,  the  target  is  fixed  at 
constant  speed  and  heading.  The  random  location  and  orientation  of  the  target  is  factored  out  of 
the  problem  by  aligning  the  coordinate  system  of  the  problem  so  that  the  target  track  lies  along 
the  x-axis.  This  is  performed  without  any  loss  of  generality,  since  the  sensor  distributions  of 
interest  are  area  coverage  systems.  (Extensions  to  barrier  problems  will  be  examined  in  a  sequel.) 
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2.1  INDIVIDUAL  SENSOR  MODEL 


For  the  high  signal-to-noise  ratio  (SNR)  application  of  interest,  each  sensor  is  essentially  an 
energy  detector  with  some  simple  classification  logic.  This  implies  that  the  detection  performance 
of  an  individual  sensor  is  dependent  only  on  the  target’s  range  to  the  sensor,  so  that  the  probability 
of  a  detection  event  for  sensor  j  is  written  as  Pd,j  =  Pd(rj),  where  Pd{r)  is  the  detection 
probability  of  a  characteristic  sensor  as  a  function  of  range  r,  and  Tj  is  the  target’s  range  to  the  jth 
sensor. 


To  examine  some  forms  of  sensor  detection  probability  functions,  several  specific  scenarios 
are  considered.  For  a  target  moving  at  a  constant  speed  vt  along  a  constant  heading,  the  path 
followed  by  the  target  over  a  period  of  time  Top  is  given  by  a  line  segment  of  length  L  =  vrTop. 
For  a  detection  to  occur,  one  scheme  may  require  that  the  sensor  lie  close  enough  to  this  line 
segment  so  that  the  target’s  closest  point  of  approach  (CPA)  to  the  sensor  r0  is  within  the 
detection  radius.  In  this  scheme,  a  detection  occurs  with  probability  Pdj  =  Pdix o)-  If  one  further 
restricts  the  problem  so  that  the  target  must  stay  above  the  detection  threshold  for  some  specified 
integration  time  Tint,  then  the  CPA  range  for  detection  opportunities  is  limited  to  (see  figure  1  for 
a  geometrical  sketch) 

to  <  yjR\  ~  {vTTint/2)\  (1) 

where  Rd  is  the  range  at  which  the  detection  threshold  is  exceeded.  Thus,  the  restriction  that 
detection  opportunities  occur  only  when  the  target  is  within  range  (above  threshold)  for  the  entire 
integration  time  yields 

f  P^i),  r,  <  /rwS  (2) 

[  0,  otherwise 


as  the  probability  of  a  single  sensor  detection  for  a  target  that  has  CPA  range  rj  to  the  sensor  j. 
This  expression  is  readily  evaluated  given  a  range-based  detection  probability,  which  may  be 
obtained  from  modeling  or  experiments.  In  either  case,  the  resulting  detection  logic  will  greatly 
impact  the  individual  sensor  performance.  For  simplicity,  one  may  wish  to  consider  a  sensor 
whose  detection  performance  versus  range  is  given  by 


Pd(rj) 


Pd,  Tj  <  Rd 
0,  Tj  >  Rd 


(3) 


Of  course,  such  an  idealized  sensor  performance  cannot  be  exactly  obtained.  To  assess  the 
performance  impact  of  deviations  from  the  ideal,  the  following  model  is  developed: 


Pd(rj;  a) 


Pd, 

(i  +  (Pd  —  Tj)/0)  Pd, 

0, 


Tj  <  Rd~  a/2 

Rd  -  a/2  <  rj  <  Rd  +  , 

Tj  >  Rd  +  a/2 


(4) 
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Figure  1.  Sensor  Detection  Geometry  Showing  Kinematic  Constraints 


where  a  >  0  is  a  parameter  to  measure  the  deviation  from  perfect,  and  the  limit  of  a  — >  0 
degenerates  to  the  case  of  the  ideal  sensor  in  equation  (3). 


2.2  FIELD  DETECTION  MODEL 

Assume  that  the  sensors  are  independently  and  identically  distributed  over  a  nominal  area 
according  to  some  known  spatial  distribution  p(x,  y).  Under  that  assumption,  the  probability  of  a 
single  sensor  having  a  CPA  range-to-target  of  less  than  r0  is  given  by 

Pr {rj  <  r0)  =  JJDP(x,y)  dxdy,  (5) 

where  D  is  the  region  of  space  given  by  D  =  {(x,  y)  :  (x  —  xt )2  +  (y  —  yr)2  <  rf}  f°r  a  target 
located  at  ( x ,  y)  =  (x^,  yr)-  For  the  case  of  a  stationary  target,  the  point  (xT,  yr)  remains 
constant  throughout  the  time  interval  of  operation  Top.  In  this  case,  the  integral  in  equation  (5)  is 
taken  as  the  area  integral  of  p(x,  y)  over  a  disk  of  radius  tq  centered  about  the  target  location. 
Specifically,  the  integral  is  given  by 

Pr (rj<r0)  =  JJDP{x,y)  dx  dy 

rro  r 2n 

—  p(r  cos(0)  +  Xt,  r  sm(6)  +  yr)  r  d6  dr,  (6) 

Jo  Jo 

where  the  transformation  to  polar  coordinates  (r,  9)  was  made  to  simplify  the  integral. 

For  the  case  of  a  moving  target  of  constant  heading  and  constant  velocity,  the  coordinate 
system  is  re-oriented  (without  any  loss  of  generality)  so  that  the  target  track  lies  along  a  line  of 
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Figure  2.  Area  D  for  Integration  Around  the  Target  Track 


constant  yr-  Then  the  integral  in  equation  (5)  becomes 


Pr(r,  <  r0)  =  JJd  p(x,  y)  dx  dy 


rt/r+To  rxT(0)-fvTTop 

/  /  p(x,  y)  dx  dy 

Jyr-ro  Jxt(  0) 

/to  rv^Top 

/  p(x'  +  xT(0 ),  y  +  yT )  dx  dy 

-ro  JO 

rrQ  p 

/  /  {p(x' +  xT{0),yT  +  r)+p{x' +  xT(0),yT-r)}  dx'dr,  (7) 

Jo  Jo 


where  the  approximation  is  valid  for  vrTop  >  ro .  The  location  of  the  area  D  in  target  track  space 
is  shown  in  figure  2.  From  the  figure,  it  is  clear  that  the  approximation  for  large  vTTop  is  to 
approximate  the  region  around  the  track  by  a  rectangle,  as  done  in  equation  (7).  In  the  case  where 
the  distance  traveled  by  the  target  ( vfTov )  is  smaller  than  the  detection  radius  ro,  a  better 
approximation  is  to  use  the  stationary  target  formula,  which  reduces  the  region  shown  in  figure  2 
to  a  circle.  A  more  exact  formula  that  holds  for  both  of  these  cases  is  readily  obtained  by 
integrating  the  exact  area  of  the  region  in  figure  2;  however,  its  evaluation  is  analytically 
intractable  and  therefore  not  included  in  this  analysis. 
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Both  of  the  formulas  (see  equations  (6)  and  (7))  obtained  for  probability  of  a  target  coming 
within  a  certain  range  tq  of  a  sensor  take  the  form 

pro 

Pr  (rj<r0)=  f{r)dr,  (8) 

J  0 

so  that  the  resulting  integrand  f(r)  in  (8)  is  a  probability  density  function  (PDF)  for  the 
range-to-target.  As  a  simple  summary  of  these  results,  recall  that  these  PDFs  are  given  by  the 
following: 

Stationary  Target  Range  PDF 

r2n 

f(r)=  /  p(r  cos(0)  +  Xt,  rsin(0)  +  yr)  rdQ,  (9) 

Jo 

Moving  Target  Range  PDF 

f(r)  =  [  T  °P  {p(x'  +  xT(0),yT  +  r)  +  p(x'  +  xT(0),yT  ~  r)}  dx',  (10) 

Jo 

where  the  target  position  is  given  by  (x,  y)  =  (xT,  yr),  which  is  either  time-varying  (second  case) 
or  fixed  (first  case). 

The  probability  of  obtaining  a  certain  range-to-target  is  combined  with  the  probability  of  the 
single  sensor  detection  at  a  given  range  to  obtain  the  probability  of  an  individual  sensor  detecting 
the  target.  The  notation  Pd<s  is  used  to  differentiate  the  detection  performance  of  an  individual 
sensor  in  the  field  (i.e.,  under  a  distribution  and  with  a  target  present)  from  the  sensor 
characteristics  that  are  described  by  as  a  function  of  range.  Since  the  detection  performance 
of  a  given  sensor  is  a  function  of  range,  and  the  range  is  a  random  variable  with  a  known  PDF,  the 
net  probability  of  detection  for  a  sensor  in  the  field  is  given  by 

Pd,s  =  /  Pd,j-f{r)dr ,  (11) 

J  0 

where  1  op  is  the  time  of  the  operational  interval  over  which  a  detection  is  sought,  and  the  infinite 
limit  will  be  reduced  by  the  sensor  performance  (i.e.,  the  sensor  has  zero  detection  capability 
beyond  a  certain  range). 

To  determine  the  net  detection  performance  of  a  field  of  sensors,  the  fact  that  the  sensors  are 
independently  and  identically  distributed  is  used.  Under  that  condition,  the  process  of  M  sensors 
out  of  N  total  sensors  detecting  the  target  is  viewed  as  a  Poisson  point  process  over  the  detection 
performance  of  an  individual  sensor.  Thus,  the  probability  of  field  level  detection  by  M  out  of  N 
sensors  is  given  by 

Pr(M  detect)  =  \  4>  =  ^Pd,s-  (12) 
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To  have  at  least  one  detection  over  the  interval  of  time  given  by  Top,  it  is  necessary  to  obtain  a 
field  level  detection  performance  of 

Pr(>  1  detect)  =  1  -  Pr(0  detect)  =  1  -  e~ <t>.  (13) 

In  a  similar  fashion,  it  is  easy  to  show  that  the  probability  of  obtaining  at  least  M  detections  over 
the  operation  interval  is  given  by 

Pr(>  M  detect)  =  1  — 


^  (<T_ 


o-4> 


(14) 


2.2.1  Special  Case  1:  Uniform  Distribution  of  Idealized  Sensors 


As  a  first  special  case,  consider  the  detection  performance  of  a  uniformly  distributed  field  of 
ideal  sensors.  In  this  case,  the  characteristics  of  an  ideal  sensor  are  that  it  has  constant  P, 'd  for 
ranges  less  than  some  range  Rd  and  that  the  Pd  is  zero  for  all  ranges  greater  than  Rd-  The 
mathematical  description  of  such  a  sensor  was  given  in  equation  (3).  Also,  the  sensor  is 
constrained  to  maintain  a  target  within  range  for  an  integration  interval  of  Pm±  as  given  in 
equation  (2).  Thus,  the  sensor  performance  is  given  by 


r<^R\-  {vrTm/W2 
r>y/R2d-  (ur7W2)2 


(15) 


The  spatial  distribution  for  sensor  location  for  a  uniform  field  over  an  area  A0  can  be  written  as 


p(x,y)  = 


V  x,y. 


(16) 


This  greatly  simplifies  the  range  PDFs  for  both  the  stationary  target  and  moving  target.  In  this 
case,  the  stationary  target  PDF  (equation  (9))  becomes 


(17) 


and  the  moving  target  PDF  (equation  (10))  becomes 


*>-jT( 


1 

Aq 


dxf  = 


2  vtT0P 
Aq 


(18) 


Combining  these  expressions  into  equation  (11)  for  the  probability  of  an  individual  sensor  in  the 
field,  and  recalling  the  identity  <f>  =  NPd,s  from  equation  (12),  the  net  sensor  field  performance  as 
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characterized  by  (f>  is  given  by 


,  2nN  fRd  , 
<p  =  — —  /  Pd-r  dr 
Ao  Jo 

Ao 


for  the  stationary  target  case,  and 

,  2 VTTopN  fy/RhlvrTintW 

< i>  =  - /  Pd  dr 

A0  Jo 


PdNQ^pUf}^,  ~  (vrTira)2 

Aq 


(19) 


(20) 


for  the  moving  target  case.  Note  that  the  cutoff  of  Rd  to  include  integration  time  was  not 
performed  for  the  stationary  case  since  it  is  a  kinematic-based  requirement,  and  the  use  of  vt  =  0 
in  the  kinematic  bound  yields  a  range  of  /?,/. 


2.2.2  Special  Case  2:  Uniform  Distribution  of  Non-Idealized  Sensors 

As  a  second  special  case,  consider  the  effects  of  a  non-idealized  sensor  detection 
characteristic,  as  shown  in  figure  3.  In  this  case,  as  described  in  equation  (4),  the  performance  of 
the  sensor  is  degraded  by  tapering  the  sensor’s  effectiveness  for  longer  ranges.  The  particular 
form  of  the  taper  is  chosen  so  that  the  total  area  under  the  detection  versus  range  curve  remains 
constant.  This  constraint  is  a  form  of  a  power  conservation  law:  the  “total”  capability  of  the 


Figure  3.  Sensor  Detection  Performance  for  Idealized  Sensor  (Dashed 
Line)  and  Non-idealized  Sensor  (Solid  Line) 
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sensor,  in  an  integrated  sense,  is  unchanged.  Also,  recall  from  equation  (4)  that  the  parameter  a 
controls  the  deviation  from  the  idealized  case,  and  a  =  0  reduces  to  the  idealized  case. 


To  examine  the  effect  of  such  a  change  in  sensor  performance,  the  analysis  in  the  previous 
section  is  repeated  for  the  new  sensor  performance  curve.  Since  the  sensor  spatial  distribution  is 
uniform,  the  stationary  and  moving  target  PDFs  are  the  same  as  in  equations  (17)  and  (18), 
respectively.  However,  the  Pd>j  term  in  the  sensor-in-field  detection  performance  expression  (11) 
has  changed  to  that  given  by  equation  (4).  Applying  this  change  to  the  sensor  field  performance 
characteristic  <j>,  one  obtains 


(j)  = 


2%N 


L 


Rd-oc/ 2  [Rd+al 2  (\  (Rd  —  r)  \ 

Pd-rdr  +  Pd  \  -  +  - - -  I  rdr 

Jr 


2 nN  Pd{Rd  -  ft/2)' 

“AT  2 


P*N*  (R2,<*2 

-A^\Ri  +  T2 


for  the  stationary  case,  and 

2 vtTopN  \ 


+  Pc 


'Rd-a/ 2  \2 

/ Rd®  «2\ 

2  12  j 


a 


=  r^pJl 

A0  [Jo  JRd-a/2  \2  a  J 


A 

2vTTopN 

A 


dr 


Pd  (■ Rd  -  | )  +  Pd  (| 


(21) 


2P(fA~ArT'0p).Rrf 

A 


(22) 


for  the  moving  target  case.  The  kinematic  requirement  that  the  target  remain  in  some  detection 
radius  for  a  given  amount  of  time  Tjnt  is  added  to  the  latter  expression  by  translating 
Rd  -4  \jR\-  (vtTm/ 2)2,  which  creates 

PdN {vTTop)  \J AR2d  -  (vTTinty 

*  = - A -  ° 


for  the  moving  target  case. 

The  moving  target  case  has  a  field-level  performance  characteristic  (j)  (see  equation  (23))  that 
is  completely  independent  of  the  taper  a.  This  occurs  whenever  the  total  area  under  the  Pd  versus 
r  curve  of  the  sensor  characteristic  remains  constant.  Also  note  that  the  stationary  target  field 
performance  characteristic  as  given  in  equation  (21)  is  always  greater  than  that  for  the  ideal 
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Figure  4.  Field  Detection  Performance  as  a  Function  of  <\> 


sensor  (a  =  0).  Since  greater  values  of  (f>  translate  to  better  detection  performance,  a  taper  on  the 
sensor  detection  characteristic  (as  long  as  the  area  is  the  same)  will  improve  performance  of  the 
system,  and  the  greater  the  taper,  the  better  the  system  performance.  Furthermore,  this 
improvement  occurs  only  for  stationary  (or  very  slowly  moving)  targets,  and  the  performance  for 
moving  targets  is  unchanged. 

For  reference,  figure  4  shows  a  plot  of  field-level  performance  characteristic  values  versus 
field-level  probability  of  detection  for  a  few  detection  limits  (i.e.,  at  least  one  detects,  at  least  two 
detect,  etc.).  For  simplicity,  the  horizontal  axis  has  been  normalized  to  be  equivalent  to 
ifra  =  N  •  Pd  for  a  uniform  distribution  over  a  100  nmi  x  100  nmi  area  with  100  m  detection 
radius  sensors.  This  provides  an  easy  reference  for  changes  in  <fr.  For  instance,  if  the  nominal  case 
has  5,000,000  sensors  with  a  sensor  Pd  of  0.8,  then  a  four-detection  requirement  yields  a  field 
detection  40%  of  the  time,  as  shown  in  the  figure  by  an  A'.  If  <f)  increases  by  25%  over  the 
nominal  case,  then  the  performance  increases  to  55%  as  shown  by  the  Y  on  the  figure. 
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2.2.3  Special  Case  3:  Non-Uniform  Distribution  of  Idealized  Sensors 


A  third  special  case  for  examination  is  the  effect  of  non-uniform  distributions  of  sensors.  To 
simplify  the  analysis,  the  discussion  is  limited  to  bivariate  Gaussian  distributions  of  idealized 
sensors.  This  case  gives  some  guidance  on  the  effects  of  sensor  clustering  (due  to  tides,  currents, 
etc.)  that  will  occur  over  time,  as  well  as  the  effect  of  rapid  deployment,  which  may  create  a 
clustered  distribution  at  the  initial  time.  As  in  the  previous  case,  the  impact  of  the  change  in 
distribution  on  the  field-level  performance  characteristic  <j>  is  examined. 


To  begin  the  non-uniform  distribution,  the  bivariate  Gaussian  distribution  is  formed  with 
independent  variables  x  and  y  as 


p(z,  V ) 


1 

2-k  o2Pq 


exp 


2  21 
-  r 

2a2 


(24) 


where  the  mean  of  the  distribution  is  taken  as  (x,  y)  =  (0,  0)  and  both  variables  are  distributed 
with  the  same  standard  deviation  o.  The  choice  to  center  the  distribution  at  the  origin  is  for 
convenience.  The  parameter  P0  is  a  scaling  term  that  represents  the  finite  extent  of  the  field,  and  it 
is  given  by 

p°  =  Hi7i^  <25) 


for  a  rectangular  area  of  dimension  L  x  L,  where  erf(-.)  is  the  error  function.  An  example  of  the 
relative  density  according  to  this  distribution  is  shown  in  figure  5.  Note  that  in  the  limit  as 
a  — >  oo,  p(x,  y)  — >  1/Ao,  which  is  the  uniform  case.  (The  asymptotic  form  of  the  error  function 
for  small  argument  is  erf(e)  «  2 e/ y/P.) 


Applying  the  bivariate  Gaussian  distribution  to  the  stationary  target  range  PDF  of  equation 
(9),  one  obtains  the  following: 

-(rcos(0)  +xt )2  -  (rsin(0)  +  j/r)2 

2^ 


f(r) 


2lU72Pc 


r2n 

r  /  exP 

0  P 


d9 


(  r  \ 

- 1 

1 

-i 

to 

1 

Is' 

^to 

+ 

^to 

< 

3 

rxTcos(0)  -  ryTsm(0) 

\2tto2pJ  6XP 

[  2cr2  J  Jo  C ^ 

(7 2 

dO 


(^)expG&) exp 


r~o4  ±vr) 

2o2 


(26) 


where  the  approximation  in  the  third  line  is  made  by  assuming  the  integrand  in  the  second  line  to 
be  nearly  equal  to  one,  which  holds  for  sensor  detection  ranges  r  that  are  small  compared  to  the 
field  distribution  standard  deviation  cr.  Applying  this  value  of  the  stationary  target  range  PDF  to 
the  sensor  field  performance  characteristic  cf),  one  obtains 


Extent  in  dimension  1  (nmi) 


Figure  5.  Sample  Spatial  Distribution  for  Bivariate  Gaussian  with  a  =  50  nmi 


By  comparing  the  value  of  the  sensor  field  performance  characteristic  4>  in  equation  (27)  with 
that  for  the  uniform  field  shown  in  (19),  one  can  assess  the  impact  of  the  non-uniformity  of  the 
field  on  the  detection  of  a  stationary  target,  which  is  parameterized  by  the  distribution  standard 
deviation  o.  The  results  of  this  comparison  are  shown  in  figure  6  for  a  variety  of  values  of  o. 
Obviously,  the  performance  approaches  the  uniform  as  cr  gets  larger  (in  the  limit  as  o  — >  oo,  one 
gets  the  same  value  for  equation  (27)  as  for  equation  (19)).  Also,  the  performance  is  better  than 
uniform  for  targets  near  the  center  of  the  distribution  (high  density)  and  worse  than  uniform  for 
those  near  the  edges  (low  density).  Of  interest  is  the  transition  from  better  than  uniform  to  worse 
than  uniform,  which  seems  to  occur  between  35  and  40  nmi  from  the  distribution  center, 
regardless  of  the  distribution  standard  deviation  a. 
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Stationary  Target 


Target  Distance  From  Distribution  Center  (nmi) 


Figure  6 .  Relative  Change  in  Field  Performance  for  a  Non-Uniform 
Sensor  Distribution  Against  a  Stationary  Target 


By  performing  a  similar  analysis  of  applying  the  bivariate  Gaussian  distribution  to  the  moving 
target  range  PDF  of  equation  (10),  one  obtains  the  following: 


f(r)  = 


2-ira2  Pa 

rvrTnp 


rVTTnp 

'-(x1  +  XT(0))2  -  (tjT  +  r)2' 

/  ^XP 
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\/2naP0 
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—  (j/r  +  r)2 
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+  exp 


-{Vt  ~ 
2  a2 


erf 


/  vTTop\ 

\2j2a)' 


(28) 


where  the  simplification  in  the  last  line  is  made  by  assuming  the  target  path  centers  about  x  =  0 
(but  makes  no  assumption  about  the  constant  yr  component  of  the  target  path).  Applying  this 
value  of  the  moving  target  range  PDF  to  the  sensor  field  performance  characteristic  4>,  one  obtains 

(using  R  =  \jB2d  -  (v-rTint/ 2)2) 


(t>  =  [  Pd-f(r)dr 

Jo 


where  P0  is  given  as  in  equation  (25). 

By  comparing  the  value  of  the  sensor  field  performance  characteristic  <fr  in  equation  (29)  with 
that  for  the  uniform  field  shown  in  (20),  one  can  assess  the  impact  of  the  non-uniformity  of  the 
field  on  the  detection  of  a  moving  target,  which  is  parameterized  by  the  distribution  standard 
deviation  a.  The  results  of  this  comparison  are  shown  in  figure  7  for  a  variety  of  values  of  a,  and 
a  target  speed  of  5  knots  over  an  observation  time  of  12  hours.  In  figure  8,  the  result  is  shown  for 
increasing  the  target  speed  to  15  knots.  One  observation  is  that  the  slow  (5-knot)  target  has 
performance  nearly  identical  to  stationary,  while  the  fast  (15-knot)  target  is  very  different.  Note 
that  for  the  fast-moving  target,  almost  all  distributions  show  performance  worse  than  uniform, 
which  occurs  because  the  target  is  moving  fast  enough  to  spend  some  of  its  time  in  the  region  of 
low  sensor  density.  The  strong  dependence  of  the  distribution  on  target  speed  is  driven  by  the 
term  vTTop  in  equation  (29).  Thus,  a  large  increase  in  observation  time  will  have  a  similar 
effect  on  the  performance  as  the  increase  in  target  speed.  Also,  the  effect  on  increased  speed  may 
be  mitigated  by  reducing  the  observation  time  appropriately. 


2.3  FIELD  FALSE  ALARM  MODEL 

For  the  field  false  alarm  analysis,  first  consider  the  characteristics  of  an  individual  sensor. 
Assume  that  the  sensors  are  independent  and  that  their  statistical  properties  of  false  alarm  are 
identically  distributed.  The  sensor  false  alarm  rate  Pfa  is  the  probability  of  an  individual  sensor 
classifying  a  detection  over  an  integration  time  interval  Tjnt  when  no  target  is  present;  that  is, 

Pfa  =  Pr {detection  during  Tjnt  |  no  target  is  present}.  (30) 


13 


Target  Moving  at  5  knots 


Figure  7.  Relative  Change  in  Field  Performance  for  a  Non-Uniform  Sensor 
Distribution  Against  a  Slow  Moving  Target 


Target  Moving  at  15  knots 


Figure  8.  Relative  Change  in  Field  Performance  for  a  Non-Uniform  Sensor 
Distribution  Against  a  Fast  Moving  Target 


The  probability  of  an  individual  sensor  registering  a  false  alarm  during  the  time  interval  of 
operation  Top  is  given  by 


Pfa,oP  =  1  -  Pr{no  false  alarm  during  Top} 
=  1-(1  -  Pfa)T°p/Tint , 


(31) 


where  the  term  (1  -  Pfa)  represents  the  probability  of  a  false  alarm  not  occurring  during  the  time 
interval  Tint. 


The  notion  of  a  field-level  false  alarm  is  not  merely  to  have  a  false  alarm  in  the  field;  rather,  it 
is  to  have  a  sequence  of  false  alarms  that  mimic  the  behavior  of  a  target.  That  is,  if  the  detection 
requirement  is  that  M  out  of  N  sensors  must  detect  within  a  region  of  interest  dictated  by 
expected  target  characteristics,  then  the  field-level  false  alarm  requirement  is  that  M  out  of  N 
sensors  must  provide  false  alarms  within  a  region  that  is  consistent  with  expected  target 
characteristics.  At  this  level  of  analysis,  one  does  not  consider  the  situation  where  the  M  sensor 
detections  are  a  mixture  of  “real”  detections  and  false  alarms.  To  account  for  the  field-level  false 
alarm  requirement,  one  forms  the  probability  of  M  false  alarms  occurring  within  a  given  region 
over  the  time  of  operation  Top. 


Given  a  sensor  spatial  distribution  of  p(x,  y ),  the  probability  of  a  single  sensor  lying  within 
the  area  D  is 

Pr(sensor  in  D)  =  Ps,d  =  JJDP(x,y)dxdy,  (32) 

where  the  entire  region  is  an  area  of  size  A0,  so  that 

JJa  p(x,  y)  dxdy  =  1.  (33) 


For  uniform  distributions,  the  probability  reduces  to  Ps,d  =  Ad/Aq,  where  Ad  is  the  area  of  D, 
which  follows  from  p(x,  y)  =  1/A0  for  the  uniform  distribution.  The  probability  of  having 
exactly  M  false  alarms  within  a  given  region  D  is  defined  by  sampling  a  Poisson  point  process 
within  region  D,  with  the  result  of  (similar  to  equation  (12)) 

(qM\ 

—  j  e~d,  0  =  NDPfa,op,  (34) 


where  No  is  the  number  of  sensors  in  region  D.  The  value  of  No  is  given  by 


Nd  —  Ps,D 


N 

Ps.Ao 


=  N  ■  Ad/ An 


(35) 


for  uniformly  distributed  fields. 
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The  problem  of  interest  for  field-level  false  alarm  analysis  is  stated  as:  What  is  the  probability 
of  at  least  M  false  alarms  occurring  within  any  area  of  size  D1  To  answer  this,  consider  the 
complementary  question:  What  is  the  probability  that  there  are  no  areas  of  size  D  that  contain  at 
least  M  false  alarms?  The  latter  question  is  more  readily  assessed  analytically.  Thus,  that 
probability  is  calculated,  and  the  complementarity  property  is  used  to  get  the  other  probability. 

The  probability  that  there  are  no  areas  of  size  D  that  contain  at  least  M  false  alarms  is  the 
same  as  asking  for  the  probability  that  every  area  of  size  D  contains  less  than  M  false  alarms. 
Beginning  with  a  single  region  of  size  D,  from  equation  (34)  one  concludes  that 

Pr(<  M  false  alarms  in 


M- 1 


- 0 


m— 0 


ml 


0  —  NpP fa, op’ 


(36) 


To  extend  this  reasoning  to  every  region  of  size  D  containing  less  than  M  false  alarms,  a  density 
argument  is  employed;  the  relative  density  of  regions  of  size  D  over  the  area  A0  is  given  by 
Ad/Aq,  so  the  number  of  them  is  A0/AD,  and  since  each  region  has  an  identical  probability  of 
less  than  M  false  alarms  (given  by  equation  (36)),  one  concludes  that 


Pr(<  M  false  alarms  in  all  regions  D )  = 


(cyNPj^)1 

m— 0 


ml 


i/a 


p  N  Pf  a, op 


(37) 


for  uniform  distributions,  where  a  =  Ad/Aq. 

Now  the  complementarity  property  is  used  to  obtain  the  final  false  alarm  expression: 


Pr(>  M  false  alarms)  =  1  — 


M- 1 

E 

Lm=0 


(ADNPfa>op/A0y 


ml 


Aq/Ad 


g  fa, op 


(38) 


which  is  the  probability  that  there  is  at  least  one  region  of  size  D  that  contains  at  least  M  false 
alarms.  The  size  of  area  D  is  dependent  on  the  detection  strategy  employed.  For  the  target 
detection  situations  analyzed  in  section  2.2,  detection  area  D  can  be  mapped  out  as  the  area  in 
which  a  sensor  has  a  detection  opportunity,  or  more  specifically 


Ad 


4>Ao 

PdN' 


(39) 


where  (f>  depends  on  whether  the  target  is  stationary  or  moving.  If  searching  for  stationary  (or 
slowly  moving)  targets,  this  yields 

AD=nRl  (40) 


and  for  faster  moving  targets, 

Ad  =  {vTTop)\fiRl  -  (vTTinty.  (41) 
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2.3.1  Special  Case  1:  Non-Uniform  Distribution  of  Sensors 


To  consider  the  effects  of  the  spatial  distribution  of  sensors  in  the  field,  a  special  case  is 
formed  where  the  field  is  broken  into  two  sections,  one  of  high  density  and  one  of  low  density. 
For  notational  simplicity,  assume  that  the  region  is  a  rectangle  in  the  x-y  plane  extending  from 
-x0  <  x  <  x0  and  -y0  <y<  yo,  for  a  total  area  of  A0  =  4 x0y0.  For  the  case  of  a  uniform 
distribution,  the  results  of  the  previous  section  can  be  applied  with  the  use  of  A o  =  4a;oPo  and 
p{x,y)  =  1/A0  for  all  (x,y). 


It  is  now  assumed  that  the  total  region  is  to  be  sectioned  into  two  regions  ( A  i  and  A2)  of  equal 
size,  where 


p(z,y) 


Pi,  -%o  <  x  <  0  (region  A\) 

2/A0  -  pi,  0  <  x  <  x0  (region  A2) 


(42) 


Since  the  area  of  each  of  the  regions  is  equal  to  2xqPq  =  Aq/2,  one  has  the  same  total  number  of 
elements  in  the  same  total  area  as  the  uniform  density  case.  For  this  special  case,  the  probability 
of  false  alarm  is  formed  as  the  probability  of  a  false  alarm  occurring  in  either  of  the  two  regions; 
to  simplify  the  model,  the  cases  of  a  false  alarm  occurring  with  some  sensors  from  each  of  the 
two  regions  are  ignored. 


The  false  alarm  probability  in  region  A\  is  given  by  equation  (38)  with  the  total  region  area 
replaced  with  the  area  of  A\  (which  is  A0/2)  and  the  total  number  of  sensors  N  replaced  with  the 
expected  number  of  sensors  in  region  A\,  which  is  given  by 


Ni  =  p{x,y)dxdy  (^—^J  =  NpiA0/2. 


(43) 


Similarly,  the  number  of  sensors  found  in  region  A  2  (which  has  the  same  area  as  A\)  is  given  by 


N2  —  f f  p{x,y)  dxdy  ( -A—\  =  N  -  NpiA0/2,  (44) 

J  J  A2  \*s,Aq/ 

where  Nx  +  N2  =  N,  as  expected.  The  probability  of  a  false  alarm  occurring  in  either  region  Ax 
or  region  A2  is  now  given  by  one  minus  the  probability  that  neither  contains  a  false  alarm;  that  is, 

Pr(>  M  false  alarms)  =  (45) 


/M- 1  M-\ 

1-  £  £ 

\m= 0  n=0 


(ZApNiPfg^p/Ao)171  {2ADN2Pfg^/A0) 


n\  ^o/(2 Ad) 


m\n 


g— NPfa,op 


The  false  alarm  performance  of  a  sample  area  distribution  relative  to  variations  of  the  relative 
size  of  Ni  and  N2  is  shown  in  figure  9.  In  this  figure,  the  horizontal  axis  represents  the  relative 
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density  of  elements  in  region  1,  where  a  relative  density  of  1  corresponds  to  the  density  for  a 
uniform  distribution  (N\  =  N2  —  N/2)  and  a  relative  density  of  2  corresponds  to  all  of  the 
elements  being  in  region  1  (N]  =  N,  N2  =  0).  The  dashed  lines  in  the  figure  correspond  to  the 
values  of  false  alarm  for  a  uniform  distribution.  For  all  three  cases  considered,  it  is  apparent  that 
small  changes  from  uniform  have  negligible  effect  on  field-level  false  alarm  performance.  Table  1 
summarizes  the  increase  in  false  alarm  rate  for  a  few  key  conditions. 


Figure  9 .  Field  False  Alarm  Rates  for  Uniform  and  Non-Uniform  Sensor  Distributions 


Table  L  Change  in  Field  Level  False  Alarm  Rate  for  Deviations  from  Uniformity 


Case 

/>,  =  1.10  x  p0 

P\  =  1.20  x  p0 

p\  =  2.00  x  pa 

N  >3 

Pfa  =  T02  X  P fa,nom 

Pfa  “  1*10  X  Pfa.nom 

Pfa  —  2.58  X  P fa,nom 

N  >4 

Pfa  ~  T06  ^  P fa,nom 

Pfa  ~  1*24  X  Pfa^nom 

Pfa~  7.64  X  P fa,nom 

N  >  5 

Pfa~  T10  X  Pfa,nom 

Pfa  =  1*40  X  Pfa,nom 

Pfa  —  15.43  X  Pfa,nom 

18 


From  table  1  it  is  clear  that  placing  all  of  the  elements  in  one  half  of  the  area  produces  large 
deviations  in  the  false  alarm  rate,  and  the  amount  of  that  change  grows  with  the  number  of 
required  detections  (i.e.,  the  effect  for  -V  >  5  is  much  larger  than  for  N  >  3).  It  is  also  clear  that 
modest  changes  from  uniformity  (on  the  order  of  10%  deviation  in  density)  have  minimal  effects 
on  the  overall  false  alarm  performance.  From  this  analysis,  it  is  concluded  that  changes  in  field 
density  of  10%  or  less  can  be  treated  as  uniform  distributions  for  false  alarm  analysis. 


2.3.2  Special  Case  2:  Clustering  of  Sensors 


The  analysis  of  fields  of  clustered  sensors  (groups  of  higher  density  in  small  areas)  can  be 
treated  as  an  extension  of  the  previous  analysis  of  the  field  broken  into  two  regions.  Consider  a 
field  made  of  C  statistically  identical  and  independent  clusters,  each  of  which  contain  N/ C 
elements.  Following  the  method  of  the  previous  section,  there  are  C  +  1  total  regions:  C  clusters 
and  one  region  of  the  remaining  empty  space.  If  the  area  of  each  cluster  is  given  by  Ac  and  there 
are  N  sensors  in  the  total  field,  then  the  probability  of  false  alarm  for  an  individual  cluster  is 
given  by 


Pr(>  M  false  alarms)  =  1  — 


M- 1 

E 

Lm=0 


{AD(N/C)Pfa,op/Acy 


-i  &cIad 


ml 


e-(N/  C)PfatOPm 


(46) 


Equation  (46)  is  the  probability  of  false  alarm  for  an  individual  cluster. 


Since  the  multiple  clusters  are  assumed  to  be  independent  and  identically  distributed  (in  a 
statistical  sense),  the  total  field  false  alarm  is  derived  from  (46)  by  use  of  the  complementarity 
property.  Recall  that  the  probability  of  no  false  alarm  in  a  given  cluster  is  given  by  one  minus  the 
probability  of  a  false  alarm.  The  field  level  false  alarm  probability  is  the  probability  that  any 
cluster  has  a  false  alarm,  which  is  given  by  one  minus  the  probability  that  all  clusters  did  not 
report  a  false  alarm.  This  is  given  by 


Pr(>  M  false  alarms)  =  1  — 
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(47) 


The  expression  for  clustered-field  false  alarm  probability  (equation  (47))  is  the  same  as  that  for  a 
uniform  field  (equation  (38))  with  the  total  area  Ao  replaced  by  the  area  of  total  cluster  coverage 
CAq ■  Note  that  this  expression  assumes  the  independence  of  clusters;  that  is,  no  consideration  is 
given  to  situations  where  a  false  alarm  occurs  because  of  sensors  from  two  neighboring  clusters 
combining  to  create  a  potential  track. 
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3.  FIELD  DESIGN  CONSIDERATIONS 


From  a  field  design  viewpoint,  the  objectives  are  low  probability  of  false  alarm  and  high 
probability  of  detection  for  a  fixed  total  number  of  sensors  N  covering  a  given  area  A0. 
Reasonable  assumptions  on  deploying  large  numbers  of  sensors  show  that  a  uniform  distribution 
is  improbable.  However,  both  the  detection  and  false  alarm  analyses  show  that  only  moderate 
deviations  from  uniformity  are  acceptable  before  performance  degrades.  For  the  examples 
considered,  one  desires  spatial  standard  deviations  on  the  order  of  tens  of  nautical  miles  (see 
figure  8)  for  detection  performance  and  density  deviations  of  no  more  than  20%  (see  figure  9)  to 
minimize  impact  on  false  alarm.  This  leads  to  the  conclusion  that  a  practical  design  solution  is  to 
deploy  a  field  of  clusters  of  sensors.  Each  cluster  can  be  held  to  a  reasonably  uniform  distribution 
by  making  the  clusters  small  enough  and  spacing  them  within  the  field  to  maximize  detection 
performance  over  the  entire  area. 

Some  ideas  for  cluster  distributions  include  widely-separated  high-density  rows  of  sensors, 
checkerboard  distributions  of  high-density  “squares”  of  sensors,  and  random  distributions  of 
clusters  of  varying  sizes.  For  the  analysis  in  this  report,  consider  the  nominal  problem  of 
searching  a  40,000-nmi2  area  for  a  contact  of  speed  up  to  5  knots.  The  probability  of  detection  for 
the  clustered  field  is  the  net  probability  of  detection  for  each  of  the  clusters;  that  is, 

Pd,field(Tc)  =  1  —  (1  -  Pd, cluster)0 1  (48) 

where  C  is  the  number  of  clusters.  The  probability  of  detection  within  a  cluster  {Pd, cluster)  is 
given  by  equation  (14),  where  N  is  replaced  by  N/C  and  <f>  is  computed  over  the  area  of  a  cluster. 
This  expression  can  be  modified  by  the  changes  shown  in  sections  2.2.1 -2.2.3  for  modifications  to 
the  ideal  sensors  in  a  uniform  distribution  within  a  cluster.  This  probability  is  only  for  the  interval 
of  time  it  takes  the  target  to  transit  a  nominal  cluster  (7c  =  ^\j  Ac /vt)-  Since  there  are  many 
such  opportunities  within  the  entire  field,  the  total  probability  of  detection  is  given  by 

PdJieldiTop)  =  1  -  (1  -  PdJield(Tc))T°’/Tc 

=  1  -  (1  -  Pd, cluster)0  T°v/Tc  •  (49) 

The  probability  of  field-level  false  alarm  for  the  clustered  field  is  given  by  equation  (47),  which 
will  hold  as  long  as  the  sensor  density  within  a  cluster  doesn’t  vary  by  more  than  10%. 

To  show  how  these  expressions  are  used  to  assist  with  field  design,  they  are  evaluated  for  a 
specific  scenario.  The  field  size  and  target  speed  are  as  given  in  the  previous  paragraph.  The 
sensors  are  all  ideal  with  a  Pd  —  0.9  out  to  a  range  of  Rd  =  100  m.  The  integration  time  is 
assumed  to  be  40  seconds  and  the  total  operation  time  under  which  a  target  is  searched  for  is 
8  hours.  Each  sensor  has  a  false  alarm  rate  of  10-5  over  the  integration  interval.  The  field  is 
populated  with  N  =  500, 000  such  sensors.  Simple  calculation  shows  the  total  detection  area 
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coverage  to  be  NitR%  =  3, 927  nmi2,  which  is  only  10%  of  the  total  search  area,  thus  creating  a 
sparse  sensor  field. 


For  this  sparse  sensor  field,  it  is  necessary  to  require  multiple  sensor  detection  reports  before 
classifying  a  detection  as  accurate.  Consider  the  effects  of  taking  N  >  3,  N  >  4,  and  N  >  5 
sensor  detection  reports  within  a  prospective  target  track.  By  evaluating  the  field-level  detection 
and  false  alarm  statistics,  the  utility  of  varying  degrees  of  clustering  is  determined.  The 
probability  of  detection  for  the  field  at  various  levels  of  clustering  (no  cluster,  C  =  100, 

C  =  1000,  C  =  10,  000)  is  shown  in  figure  10  as  a  function  of  cluster  size.  It  is  obvious  that 
detection  performance  is  always  less  for  clustering  than  for  the  uniform  field.  Also  note  that  each 
value  of  C  (number  of  clusters)  has  a  corresponding  “sweet  spot”  in  cluster  size.  As  long  as  the 
cluster  is  sized  according  to  this  sweet  spot,  detection  performance  is  nearly  identical  for  the 
various  cluster  sizes.  However,  the  width  of  the  sweet  spot  (how  tightly  held  cluster  size  needs  to 
be)  decreases  with  larger  numbers  of  clusters.  This  makes  the  most  practical  solution  to  have  a 
relatively  smaller  number  of  clusters  and  keep  the  cluster  size  in  the  sweet  spot. 


Figure  10.  Detection  Performance  for  a  Field  of  Area  40,000  nmi2 
(Pink:  Uniform  Distribution,  Blue:  C=20,  Green:  C=100,  Red:  C=500) 


The  motivation  for  clustering  was  to  decrease  the  false  alarm  rate  for  the  field.  The  effects  of 
clustering  on  field  false  alarm  rate  are  shown  in  figure  1 1  for  the  same  levels  of  clustering  (no 
cluster,  C  =  100,  C  =  1000,  C  =  10, 000).  Here  it  is  obvious  that  clustering  improves  the  false 
alarm  performance  (compared  to  uniform)  if  the  clusters  are  made  large  enough.  Furthermore,  the 
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Cluster  Radius  [nmi] 


Figure  11.  False  Alarm  Performance  for  a  Field  of  Area  40,000  nmi 2 
(Pink:  Uniform  Distribution,  Blue:  C=20,  Green:  C=100,  Red:  C=500) 


cluster  size  must  be  larger  for  smaller  numbers  of  clusters.  Thus,  a  practical  guideline  may  be 
derived  that  the  cluster  size  be  chosen  as  the  largest  cluster  where  the  density  can  be  held 
relatively  constant  (within  10%  variation).  Using  this  guideline,  the  number  of  clusters  is  chosen 
so  that  cluster  size  lands  in  the  “sweet  spot”  for  detection  performance.  In  this  way,  a  reasonably 
optimal  tradeoff  is  made  between  false  alarm  performance  and  detection  performance.  Obviously, 
these  decisions  must  be  made  with  regard  to  meeting  any  set  absolute  performance  criteria  (i.e., 

P,i  must  be  greater  than  X,  and  P/a  must  be  less  than  Y)\  however,  this  analysis  shows  some  of 
the  tradeoffs  that  can  be  advantageously  made. 


4.  CONCLUSIONS 


This  report  has  provided  an  analytical  derivation  of  the  important  tradeoffs  involved  in  the 
performance  of  large  distributed  systems  of  proximity  sensors.  Throughout  the  analysis,  realistic 
features  were  added  to  examine  their  impact  on  modifying  the  idealized  system’s  resulting 
performance.  A  summary  of  the  important  results  follows: 
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•  Sensor  detection  performance  versus  range  is  driven  by  the  total  area  under  a  Pd 
versus  range  curve. 

•  Field  detection  performance  for  moving  targets  is  generally  better  than  for  stationary 
targets,  especially  when  an  M  out  of  N  detection  criterion  is  applied  to  determine 
prospective  targets. 

•  Field  probability  of  detection  is  lowered  by  increasing  the  number  of  required 
detections;  however,  probability  of  false  alarm  is  lowered  by  a  greater  amount. 

•  Non-uniformity  of  the  spatial  distribution  degrades  field  detection  performance  in 
general.  The  degradation  is  larger  for  faster  moving  targets  than  for  slow-moving  or 
stationary  targets.  Also,  the  degradation  is  not  very  severe  if  the  standard  deviation  of 
the  distribution  is  on  the  order  of  the  field  size. 

•  Non-uniformity  of  the  spatial  distribution  degrades  field  false  alarm  performance  in 
general.  The  degradation  is  moderate  if  the  field  sensor  density  does  not  vary  by  more 
than  10%.  The  degradation  with  non-uniformity  is  worse  for  larger  numbers  of 
required  detections. 

•  Clustering  of  sensors  provides  an  effective  way  to  lower  field  probability  of  false 
alarm  at  the  expense  of  reduced  field  probability  of  detection.  However,  the  effect  on 
detection  can  be  minimized  by  choosing  the  cluster  size  and  number  of  clusters 
carefully  to  lie  in  a  detection  “sweet  spot.” 

•  Clustering  of  sensors  provides  a  means  of  mitigating  the  effects  of  non-uniformity, 
since  the  degree  of  non-uniformity  is  important  only  within  a  cluster.  Thus,  since 
clusters  are  geometrically  small,  their  distributions  are  easier  to  control  within 
acceptable  bounds. 


These  guidelines  and  the  supporting  calculations  show  that  it  is  possible  to  optimize  the 
performance  of  a  sensor  field  for  a  given  sensor  type  by  choosing  a  proper  clustering  scheme, 
distribution  requirement,  and  number  of  detections  to  generate  a  track.  It  is  intended  that  the 
analysis  contained  in  this  report  be  used  and  expanded  to  assist  in  the  design  of  optimally 
performing  systems  of  large  numbers  of  distributed  sensors.  Extensions  that  should  be  done 
include  combined  analysis  with  tracking  systems,  analysis  of  the  effects  of  mixed  (real  plus  false 
alarm)  reports,  analysis  of  fields  of  mixed  sensor  type,  and  numerical  optimization  of  field 
performance  in  a  Pareto  optimal  sense. 
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